Subharmonic resonance of global climate to solar forcing 
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It is shown that, the wavelet regression detrended fluctuations of the monthly global temperature 
data (land and ocean combined) for the period 1880-2009yy, are completely dominated by one-third 
subharmonic resonance to annual forcing (both natural and anthropogenically induced). Role of 
the oceanic Rossby waves and the resonance contribution to the El Nino phenomenon have been 
discussed in detail. 

PACS numbers: 92.70.Gt, 92.70.Qr, 92.10.am, 92.10.Hm 
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INTRODUCTION 



The monthly global temperature data are known to 
be strongly fluctuating. Actually, the fluctuations are 
of the same order as the trend itself (see figures 1 and 
2). While the nature of the trend is widely discussed 
(in relation to the global warming) the nature of these 
strong fluctuations is still quite obscure. The problem 
has also a technical aspect: detrending is a difficult task 
for such strong fluctuations. In order to solve this prob- 
lem a wavelet regression detrending method was used in 
present investigation. Then a spectral analysis of the 
detrended data reveals rather surprising nature of the 
strong global temperature fluctuations. Namely, the de- 
trended fluctuations of the global temperature for the last 
century (the instrumental monthly data are available at 
Ref. Ml see also Ref. (2|) are completely dominated 
by so-called one-third subharmonic resonance to annual 
forcing. 
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FIG. 1: The monthly global temperature data (dashed line) 
for the period 1880-2009. The solid curve (trend) corresponds 
to a wavelet (symmlet) regression of the data. 
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FIG. 2: The wavelet regression detrended fluctuations from 
the data shown in Fig. 1. 



SUBHARMONIC RESONANCE TO ANNUAL 
FORCING 



There are many well known reasons for asymmetry in 
response of the North and South Hemispheres to solar 
forcing: dominance of water in the Southern Hemisphere 
against dominance of land in the Northern one, topo- 
graphical imbalance of land (continents) and oceans in 
the Northern Hemisphere due to continental configura- 
tion, seasonality and vegetation changes are much more 
pronounce on land than on ocean surface, and anthro- 
pogenically induced asymmetry of the last century. This 
asymmetry results in annual asymmetry of global heat 
budget and, in particular, in annual fluctuations of the 
global temperature. Nonlinear responses are expected as 
a result of this asymmetry. 

Figure 1 shows (as dashed line) the instrumental 
monthly global temperature data (land and ocean com- 
bined) for the period 1880-2009, as presented at the 
NOAA site [l[ . The solid curve (trend) in the figure cor- 
responds to a wavelet (symmlet) regression of the data 
(cfRefs. (3,0). Fi gure 2 shows corresponding detrended 
fluctuations, which produce a statistically stationary set 
of data. Most of the regression methods are linear in 
responses. At the nonlinear nonparametric wavelet re- 
gression one chooses a relatively small number of wavelet 
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FIG. 3: Spectrum of the wavelet regression detrended fluctu- 
ations of the monthly global temperature anomaly (land and 
ocean combined). 



FIG. 5: The same as in Fig. 3 but for the wavelet regression 
detrended fluctuations of the global Sea Surface Temperature 
Anomalies. 
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FIG. 4: The monthly global Sea Surface Temperature Anoma- 
lies (dashed line) for the period 1850-2008yy (the data taken 
from Ref.[l7|). The solid curve (trend) corresponds to a 
wavelet (symmlet) regression of the data. 



coefficients to represent the underlying regression func- 
tion. A threshold method is used to keep or kill the 
wavelet coefficients. In this case, in particular, the Uni- 
versal (VisuShrink) thresholding rule with a soft thresh- 
olding function was used. At the wavelet regression the 
demands to smoothness of the function being estimated 
are relaxed considerably in comparison to the traditional 
methods. Figure 3 shows a spectrum of the wavelet re- 
gression detrended data calculated using the maximum 
entropy method (because it provides an optimal spectral 
resolution even for small data sets). One can see in this 
figure a small peak corresponding to a one-year period 
and a huge well defined peak corresponding to a three- 



years period. 

In order to understand underlying physics of the very 
characteristic picture shown in the Fig. 3 let us imagine 
a forced excitable system with a large amount of loosely 
coupled degrees of freedom schematically represented by 
Duffing oscillators (which has become a classic model for 
analysis of nonlinear phenomena and can exhibit both 
deterministic and chaotic behavior depending on 

the parameters range) with a wide range of the natural 
frequencies ujq (it is well known Q that oscillations with 
a wide range of frequencies are supported by ocean and 
atmosphere, cf also Ref. [loj): 



x + cjq x + jx + (3x 3 = F sinuit 



(1) 



where x denotes the temporal derivative of x, (3 is the 
strength of nonlinearity, and F and u> are characteristic 
of a driving force. It is known (see for instance Ref. [llj ]) 
that when u> ~ 3u; and ^«1 the equation (1) has a 
resonant solution 



x(t) 



a cos 



F 



COS U)t 



(2) 



where the amplitude a and the phase ip are certain con- 
stants. This is so-called one-third subharmonic resonance 
with the driving frequency uj corresponding to the annual 
NS-asymmetry of the solar forcing (the huge peak in Fig. 
3 corresponds to the first term in the right-hand side of 
the Eq. (2)). For the considered system of the oscillators 
an effect of synchronization can take place and, as a con- 
sequence of this synchronization, the characteristic peaks 
in the spectra of partial oscillations coincide 12| . It can 
be useful to note, for the global climate modeling, that 
the odd-term subharmonic resonance is a consequence of 
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ROLE OF ROSSBY WAVES 
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FIG. 6: Defect of autocorrelation function versus r (in ln-ln 
scales) for the wavelet regression detrended fluctuations of the 
global Sea Surface Temperature Anomalies. The straight line 
is drawn in order to indicate scaling Eq. (3). 
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FIG. 7: Spectrum of sea surface height fluctuations [TH ] 
(TOPEX/ Poseidon and ERS-1/2 altimeter measurements) 
in the logarithmic scales. The profound peak corresponds to 
the annual cycle. The dashed straight line is drawn in order 
to indicate correspondence to the scaling Eq. (3). 



The fluctuations of oceanic temperature cause certain 
variations of the sea surface height. These variations are 
intermixed with the sea surface height variations caused 
by the oceanic planetary Rossby waves. The oceanic 
planetary Rossby waves play an important role in the 
response of t he g lobal ocean to the forcing (see, for in- 
stance, Refs. [lj], [lj|) and they are of fundamental im- 
portance to ocean circulation on a wide range of time 
scales (it was also suggested that the Rossby waves play 
a crucial role in the initiation and termination of the 
El Nino phenomenon, see also below). Therefore, they 
present a favorable physical background for the global 
subharmonic resonance. For that reason it is interesting 
to look also separately at global sea surface temperature 
anomalies. These data for time range 1850-2008yy are 
shown in Fig. 4 (the monthly data were taken from Ref. 
[l7j | . see also Ref. [Hj]). The solid curve (trend) in the 
figure corresponds to the wavelet (symmlet) regression 
of the data. Figure 5 shows a spectrum of the wavelet 
regression detrended data calculated using the maximum 
entropy method. The spectrum seems to be very similar 
to the spectrum presented in Fig. 3. 

Since the high frequency part of the spectrum is cor- 
rupted by strong fluctuations (the Nyquist frequency 
equals 0.5 \mon~ 1 ]), it is interesting to look at corre- 
sponding autocorrelation function C(r) in order to un- 
derstand what happens on the monthly scales. It should 
be noted that scaling of defect of the autocorrelation 
function can be related to scaling of corresponding spec- 
trum: 

l-C(r)~T a E{f)~f- {1+a) (3) 

Figure 6 shows the defect of the autocorrelation func- 
tion in ln-ln scales in order to estimate the scaling ex- 
ponent a ~ 0.6 ± 0.04 (the straight line in this fig- 
ure indicates the scaling Eq. (3)). The existence of 
oceanic Rossby waves was confirmed rather recently by 
NASA/CNES TOPEX/Poseidon satellite altimetry mea- 
surements. Corresponding to these measurements spec- 
trum of the sea surface height fluctuations, calculated in 
Ref. [HI with daily resolution (see also [16]), is shown 
in Figure 7. The dashed straight line in this figure is 
drawn in order to indicate correspondence to the scaling 
Eq. (3): 1 + a ~ 1.6. 

EL NINO PHENOMENON 



the reflection symmetry of the natural nonlinear oscilla- 
tors (inyariance to the transformation x — > — x, cf. also 
Ref. [H). 



The Rossby waves (together with Kelvin waves) and 
a strong atmosphere-ocean feedback provide physical 
background for the El Nino phenomenon (see, for 
instance, Ref. [l~9} and references therein). Figure 9 
shows spectrum for the wavelet detrended fluctuations 
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FIG. 8: The monthly Global-SST ENSO index (dashed line) 
for the period 1850-2008yy (the data taken from Ref. 
the index is in hundredths of a degree Celsius) . The solid 
curve (trend) corresponds to a wavelet (symmlet) regression 
of the data. 
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FIG. 9: The same as in Fig. 3 but for the wavelet regression 
detrended fluctuations of the Global-SST ENSO index. 



of the so-called Global-SST ENSO index (Fig. 8), 
which captures the low-frequency part of the El Niho 

jhenomenon (the monthly data are available at Ref. 

20]). The annual forcing can come from the oceanic 
Rossby waves (cf Fig. 7). To support this relationship 
we show in figure 10 defect of autocorrelation function 
calculated using the wavelet detrended fluctuations from 
Fig. 8. The ln-ln scales have been used in Fig. 10 in 
order to estimate the scaling exponent a ~ 0.6 ± 0.03 
(the straight line in this figure indicates the scaling Eq. 
(3), cf Figs. 6 and 7). Using these observations one can 
suggest that the El Niho phenomenon has the one-third 
subharmonic resonance as a background. 
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FIG. 10: The same as in Fig. 6 but for the wavelet regression 
detrended fluctuations of the Global-SST ENSO index. 



The data were provided by National Climatic Data 
Center at NOAA and by Joint Institute for the Study of 
the Atmosphere and Ocean. I also acknowledge that a 
software provided by K. Yoshioka was used at the com- 
putations. 
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